Exact SO(5) Symmetry in spin 3/2 fermionic system 
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The spin 3/2 fermion models with contact interactions have a generic S0(5) symmetry without 
any fine-tuning of parameters. Its physical consequences are discussed in both the continuum and 
lattice models. A Monte-Carlo algorithm free of the sign problem at any doping and lattice topology 
is designed when the singlet and quintet interactions satisfy Uq < U2 < —^Uo {Uo < 0), thus making 
it possible to study different competing orders with high numerical accuracy. This model can be 
accurately realized in ultra-cold atomic systems. 
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With the rapid progress in ultra-cold atomic systems, 
many alkali fermions have been cooled below Fermi tem- 
peratures 

[1 H 13' All of them except ^Li have spms 
higher than 1/2 in the lowest hyper-fine multiplets. The 
spin degrees of freedom become free in the optical traps, 
which has attracted interest in their effects on Cooper 
pair structures and collective modes 0, |^ . The proposal 
of the optical lattice 6] has led to a tremendous progress 
in studying the strongly correlated bosonic lattice sys- 
tems Recently, fermionic lattice systems 
are also exciting. For example, the degenerate ''"K gas 
has been prepared in a one-dimensional optical lattice 

m 

Comprehensive analysis of symmetries is helpful in un- 
derstanding the physics in stroiigly correlated systems. 
For example, the SO (5) theory of high Tc cuprates 
unifies the d-wave superconductivity (SC) and antiferro- 
magnetism (AF) orders, leading to many experimental 
consequences. The sharp neutron scattering mode can 
be interpreted as the pseudo-Goldstone mode[la|. The 
prediction of the antiferromagnetic vortex core [lj| has 
also been verified in recent experiments p^ . 

In this article, we focus on the symmetry properties 
and corresponding consequences in the spin 3/2 system 
with contact interactions, including both the continuum 
model with s-wave scattering and the generalized lat- 
tice Hubbard model with on-site interactions. For neu- 
tral atoms, these interactions are generally described by 
two parameters in the total spin St 0, 2 channels as 
90,2 = 4:TTh'^ao,2/M in the continuum model with 09,2 the 
corresponding s-wave scattering lengths and M the atom 
mass; or Uo,2 in the lattice model. Interactions in the odd 
total spin {St — 1, 3) channels are forbidden by Pauli's 
exclusion principle. Remarkably, in addition to the ex- 
plicit spin SU(2) symmetry, an enlarged S0(5) symmetry 
is present without any fine tuning of parameters. In the 
continuum model, this symmetry has direct consequences 
on the collective modes and pairing structures. In the lat- 
tice model, exact phase boundaries of various competing 
phases can be determined directly from symmetries. Be- 
cause of the time-reversal symmetry of the Kramers dou- 
blets, a Monte-Carlo algorithm free of the notorious sign 
problem is designed when Uq < U2 < —3/5 Uo {Uq < 0) 



at any filling level and lattice topology. 

We start with the standard form of the spin 3/2 Hamil- 
tonian of the continuum model 0, |^ 



a=±3/2,±l/2 
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M)V'a(r) 



+ 5oPj,o(r)Fo,o(r)+52 ^ 4™W^2,™(r)}, (1) 



m=±2,±l,0 



with d the space dimension, the chemical potential and 
Plo,Plm the singlet {St = 0) and quintet {St = 2) 
pairing operators defined through the Clebsh-Gordan co- 
efficient for two indistinguishable particles as Pp ^ (r) = 

E„/3(||;F,m|f|a/?)Vt(r)v4(r), where F = 0, 2 and 
m = -F,-F + 1,...,F. 

We first construct the SO (5) algebra by introducing 
the five Dirac F" (1 < a < 5) matrices 
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Then the ten SO (5) generators are defined as F"'' ~ 
-|[F",r''](l < a,b < 5), where / and a are the 2x 2 
unit and Pauli matrices. The four-component spinor can 

be defined by '0(r) — (V'a (r), V'i (r), V'-i (r), "0-^ (i")) ■ 
Furthermore, the bilinear operators can be classified ac- 
cording to their properties under the SO (5) transforma- 
tions. The 16 bilinear operators in the particle-hole (p-h) 
channel can be classified as S0(5)'s scalar, vector, and 
anti-symmetric tensors (generators) as 



n(r) = Vl(r)0a(r), n,(r) = l^t (r)r-^^^(r), 



L„fc(r)---Vl(r)PS>^(r 



(2) 



Lab and ria together form the SU(4), or isomorphically, 
the S0(6) generators. The spin SU(2) generators Jx,y,z 
are expressed a.s J+ ~ Jx + i Jy ~ \/3(— L34 + 1^24) + 

(ii2 + 1^25) - j(ii3 + iL^b), J- ^ J+, and = -L23 + 
2 ii5. n and ria have spin and 2, and Lab contains 
both the spin 1 and 3 parts. Pairing operators can also 
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be organized as S0(5) scalar and vectors through the 
matrix R = TiT^ 

ry^(r) = Rerj + i Imr] = ^'4!l{r)Rai3'4'l{r), 

Xl(r) = Rcxa + tlmxa = -^^i(r)(r'^i?)„;3V'^(r)(3) 

where Pj,o = -^VV2, and pI±^ = {tx\ + «xi)/2, 
pI±i = (-X^ ± *xD/2, 4o = -*xi/V2. That is, xl 
are polar combinations of J^'s eigenoperators Pj m- The 
existence of the R matrix is related to the pseudoreality of 
S0(5)'s spinor representation. It satisfies i?^ = —1, R^ = 
= *R = -R and RV'R = ~*r'^, RV'^R = *r'''' 
|16|. The anti- unitary time-reversal transformation can 
be expressed as T — R C, where C denotes complex 
conjugation and — —1. N, Ua, and Lab transform 
differently under the T transformation 

TnT^^ = n, TnaT~^ = Ha, TLabT^^ = -Lab- (4) 

With the above preparation, the hidden SO (5) sym- 
metry becomes manifest. The kinetic energy part has 
an explicit SU(4) symmetry which is the unitary trans- 
formation among four spin components. The singlet 
and quintet interactions are proportional to 77^(r)r/(r) 
and xl(r)xa(r) respectively, thus reducing the symme- 
try group from SU(4) to SO (5). When go = g2, the 
SU(4) symmetry is restored because xL together form 
its 6 dimensional antisymmetrical tensor representation. 
In the continuum model, interactions in other even par- 
tial wave channels also keep the SO (5) symmetry. The 
odd partial wave scattering include spin 1 and 3 channel 
interactions gi and (73, which together could form the 10- 
d adjoint representation, of SO (5), if and only gi — gs- 
However, to the leading order, p-wave scattering is much 
weaker than the s-wave one for neutral atoms, and can 
thus be safely neglected. 

The SO (5) symmetry implies more degeneracies in the 
collective excitations in the spin 3/2 Landau fermi liquid 
theory, which generally requires four fermi liquid func- 
tions in total spin St = 0,1,2,3 channels. The S0(5) 
symmetry of the microscopic Hamiltonian reduces these 
to three independent sets, classified according to the 
SO (5) scalar, vector, and tensor channels as 

UpMip^p') = fs{p,p') + Mp,p')irV2)^i,{Ty2)^s 

+ Mp,p'W'/2UiT^'>/2),s. (5) 

In other words, the effective interaction functions in the 
St = 1,3 channels are exactly identical in all orders 
in perturbation theory. Furthermore, within the s-wave 
scattering approximation, the interaction functions be- 
come constants, and are given as fs = {go + 5172)/ 16, 
fv = (ffo - 352)74, ft = -{go+ S'2)/4- Experiments in 
the fermi liquid regime can determine the four fermi liq- 
uid constants in the St = 0,1,2,3 channels separately 
and verify the degeneracy between spin 1 and 3 chan- 
nels. This degeneracy appears to be accidental in Ref. 



. It is in fact exact and protected by the generic S0(5) 
symmetry. 

The S0(5) symmetry also enriches the Cooper pair 
structures. Ref. showed that, in addition to the 
singlet pairing when go < 0, the spin 3/2 system en- 
ergetically favors the polar pairing state in the quin- 
tet channel when (72 < with the order parameter 

Keai = 6^0 + 6(^2 + Pl-2) « 6x1 + 6x1 (6,2 are 
real) . We understand that this is only a special case of the 
general pairing structures spanned by all the Xi~5- The 
polar pairing states break the SO(5)(8) U(l) (charge) sym- 
metry to 50(4) (Xi Z2, and thus the Goldstone (GS) man- 
ifold 'is the quotient space [SO{5)(g)U{l)]/[SO{'i)(g)Z2] = 
[S** (8) U{l)]/Z2. Its dimension, 5, is the number of GS 
modes. When both 50,2 are positive, s-wave pairing is 
not favorable. However, similarly to the spin fluctuation 
exchange mediated p-wave pairing in 3He [Hill, the 
spin fluctuations in the tensor channel can induce effec- 
tive attractions between two atoms with total spin 1 and 
3. This may lead to p-wave pairing where the spin part 
forms the 10-d adjoint representation of S0(5). 

Now we consider the more interesting case of spin 3/2 
fermions in the optical lattice. The periodic potential 
is V{x,y,z) = Vb(sin^(fca;) + sin^(fcy) -I- sin^(fcz)) with 
Vo the potential depth , k — tt/Zq the wavevector, and 
lo the lattice constant. The hopping integral t between 
neighboring sites decreases exponentially with increasing 
Vq. Within the harmonic approximation, the parame- 
ter U/AE w {TT^/2){as/lo){Vo/Ery^'^, with U the repul- 
sion of two fermions on one site, AE the gap between 
the lowest and first excited single particle state in one 
site, as the s-wave scattering length in the correspond- 
ing channel, and Er = fi^k'^ /2M the recoil energy. With 
the typical estimate of ~ lOOas ( as the Bohr ra- 
dius), lo - 5000A, and {Vo/Erfl'^ « 1 - 2, we arrive at 
U/AE < 0.1. Thus this system can be approximated by 
the one-band Hubbard model 

H = -tY^ {c\^Cj„ 4- h.c] -f^J2^^'^^"^ 

+uo J2 pU^)PoA^) + U2 plmi^)P2,r..{m) 

i i,m=±2,±l,0 

for particle density n < 4. At half-filling on a bipar- 
tite lattice, fi is given by ^0 = {Uo + 5[/2)/4 to en- 
sure the particle-hole (p-h) symmetry under the trans- 
formation Ci^a — > {—Ycl^. The lattice fermion oper- 
ators and their continuum counterparts are related by 
'0a(r) = Ca{i) / {loY^^ ■ We use the same symbols for bi- 
linear fermion operators as in the continuum model. 

The proof of SO (5) invariance in the continuum model 
applies equally well in the lattice model at any lattice 
topology and at any filling level. Eq. |^ can be conve- 
niently rewritten in another manifestly S0(5) invariant 
form as 



i 



3 



Hr = 



E 



j.l<a<5 



16 



- 2) n, 



(8) 



where the SU(4) symmetry appears at Uq — U2 as before. 



Phase order parameters 
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TABLE I: Order parameters, the corresponding Goldstone 
manifolds and the number of Goldstone modes in each phase 
on the bipartite lattice at half-filling. 



MC boundary 
3 Uo= -5U2 



F:S0(7)line 
Uo= -3U2 



C: < Tj(i) >^ 



G: SO(5)*SU(2) line 
Uo= 5 U2 



MC boundary 

U0=U2 



E: SU(4) line 
Uo=U2 



B:<(-)'«„(0>^0 




-''D: <(-)'#(/) >^0 



H: SO{7) line 
Uo= -3 U2 



FIG. 1: The mean-field phase diagram at half-filling on a 
bipartite lattice. A) and B): staggered phases of the SO(5) 
adjoint and vector Reps; C): the singlet superconductivity; 
D): CDW; E), F), G) and H): exact phase boundaries with 
higher symmetries. Between the dashed lines {Uo < U2 < 
—3/5 U2), a Monte-Carlo algorithm free of the minus sign 
problem at any filling level and lattice topology is possible. 



The lattice Hamiltonian, Eq. 1^1 contains even higher 
symmetries under certain conditions. One can construct 
the largest SO(8)|l3 algebra using all the independent 
fermionic bilinear operators. Its generators Mat (0 < 
a < b < 7) including Lab(l < a < 6 < 5) as its S0(5) 
sub-algebra, are denoted as 



Mab = 



/O Rexi - Rex5 N 




Imxi 


ni 


Lab ^ 




Imxs 


"5 





— Im?7 


V 


1 



with N = (n ~ 2)/2. Its Casimir is a constant Csos = 
Eo<a<b<7Mibi^) = 7. The global S0(8) generators 
are defined to be uniform in the p-h channel as Mab = 
J2i^ab{i) and staggered in the p-p channel as Mab = 
)*Mafc(i) on the bipartite lattice. These global gen- 
erators commute with the hopping term [Mab, Hq] = 0. 
On the other hand, order parameters transformed under 
the SO (8) group should be staggered (uniform) in the 
p-h (p-p) channel respectively. The SO (8) symmetry is 
always broken by the interaction, but its subgroup sym- 
metry, the S0(5)® SU(2) and S0(7) symmetries, appear 
under special conditions as shown below. 

At Uq = 5 U2, Hi can be rewritten as Hj — 
Y.^,l<a.b<b{~U2 Ll^{i) - (m - Aio)n(z)}, using the Fierz 

identity Ei<a<5ia6(*) + Ei<a<5 "a W + 5iV^'(i) = 5- As 
a generalization of the pseudospin algebra in the usual 
Hubbard model |23|, we construct them as ryt, rj, N . The 
symmetry at half-filling is S0(5) ® SU(2), which unifies 
the charge density wave (CDW) and the singlet pairing 



(SP) order parameters. Away from half-filling, this sym- 
metry is broken but 77, ry^ are still eigen-operators since 
[H, r^t] = -(/X - /xo)r?t, and [H, rj]^{^- ^o)?7- 

The p-h channel S0(5) (8) U(l) symmetry can also be 
extended to S0(7) at Uq — —3U2 where Hj can be rewrit- 
ten as Hi ^ E*,o<a<b<6 {§^^2 Mab(j)^ - (a* - Aio)"(i)}- 
The SO (7) symmetry is exact at half-filling. Its 7-d 
vector Rep unifies the staggered 5-vector and SP or- 
der parameters. Its 21-d adjoint Rep unifies the stag- 
gered SO (5) adjoint Rep order parameters, CDW, and 
quintet pairing (QP) order parameters. Away from half- 
filling, QP operators are spin 2 quasi- Goldstone operators 
W, xU = -(a* - Mo)xi and [H, Xa] = Mo)Xa- These 
X modes are just the analogs of the tt modes in the high 
Tc context 13]. 

The above symmetry structures guide the mean field 
(MF) analysis of the phase diagram. In the weak coupling 
limit, the complete MF decoupling is performed in the 
direct, exchange, and pairing channels. We take the MF 
ansatz on the 2d square lattice 

M)) - {-yna_, {N{^)) = i-yN, 

{Lab{l)) = {-yLab m)=% (Xa(*))-Xa, (9) 

and solve it self-consistently at half-filling to obtain the 
phase diagram shown in Fig. 1. Higher symmetry lines E, 
F, G, H separate phases A, B, C, D as first order phase 
transition boundaries where order parameters smoothly 
rotate from one phase to another. Symmetries on lines 
E, F, G, H and the order parameters are SU(4) (adjoint 
Rep), SO (7) (vector Rep), S0(5)(g) SU(2) (scalar ® vector 
Rep), SO (7) (adjoint Rep) as discussed before. Phases A 
and B spontaneously break the SO (5) symmetry in the 
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adjoint and vector Rep channels respectively. Phases C 
and D have singlet pairing SC and CDW as order pa- 
rameters, respectively. Order parameters in each phase 
and corresponding GS modes are summarized in Table 
1. The effective theory is generally given by a quantum 
non-linear a model defined on the GS manifold. 

One major difficulty of Monte-Carlo simulations in 
fermionic systems, the sign problem 21], is absent in the 
spin 3/2 model when Uq < U2 < —3/5 Uq. By the 
Hubbard-Stratonovich (HS) transformation, the parti- 
tion function can be written as below when V = — (3C/o + 
5U2)/8 > and = {U2 - Uo)/2 > 0, or equivalently 
Uo<U2< -3/5 Uo, 

Z = J Dn J Dn'' exp^-^ J dT^n{i,Tf 
- Y f dTY,nl{i,T)] det{/ + B}, 

^ i,a 

„ ^ „ ffo-i--H'/(i-) J- ^j^g time order 
operator. Its discrete version is 

Hi{r) - - (^(*'^) - 2) + (m - Mo)} 

■i 

n«(*,r), (10) 

where At — (3/L. I -\- B is invariant under the time- 
reversal transformation: T{I + B)T^^ = I + B. If A is 
an eigenvalue of / -f -B with the eigenvector then A* 
is also an eigenvalue with the eigenvector T \(f)). From 
T2 = -1, it follows that ((l)\T(t)) = {T^(l)\T(j)) = 0, i.e. |0) 
and T\<j)) are orthogonal. Thus although I + B may not 
be Hermitian because of the T operator, its determinant, 
a product of A* A, is always positive semi-definite. Our 
proof is equally valid in the practical sampling with the 
discrete HS transformation as in Ref. [23|, and has been 
confirmed numerically [22. We emphasize that this proof 
is valid for any filling and lattice topology. A similar 
model has recently been introduced in Ref. p^ . where 
the sign problem is also absent. However, their model 
keeps only the diagonal n1 interaction and is not spin 
rotationally invariant. The valid region for the above al- 
gorithm (see Fig. 1) includes the 5-vector phases B, SP 



phase C and their S0(7) boundary, which are analogs of 
the competitions between antiferromagnetism and super- 
conductivity in the high Tc context. It would be inter- 
esting to study the doping effect, the frustration on the 
triangular lattice, etc, which are difficult at low tempera- 
tures for previous Monte-Carlo works. Extensive numer- 
ical simulations are currently being carried outj24|. 

Besides the alkali atoms, the trapping and cooling of 
the alkaline-earth atoms are also exciting recently ^25.. .26j. 
Among these two famihes, ^'^^Cs, ^Be, ^^^Ba and ^^^Ba 
are spin 3/2 atoms. ^'^^Cs is unstable and the 2s^ 
2s^2p^ resonance of Be lies in the ultraviolet region, mak- 
ing them difficult for experimental use. The resonances 
of the last two Ba atoms are 6s^ — > 6s^6p^ at 553.7 nm 
[27l |. thus they are possible candidates. Their scattering 
lengths are not available now, but that of ^^^Ba (spin 0) 
was estimated as — 41aB Because the 6s shell of Ba 
is full-filled, both the oo, 02 of ^^^Ba and ^^"^Ba should 
have the similar value. Considering the rapid develop- 
ment in this field, we expect more and more spin 3/2 
systems will be realized experimentally, allowing us to 
explore the full phase diagram. 

In summary, we found an exact and generic SO (5) sym- 
metry in spin 3/2 models with local interactions. This 
model can be accurately realized in cold atom systems 
and the theoretical predictions can be tested experimen- 
tally by the exact relationship among the Landau fermi 
liquid parameters, spin wave functions of the Cooper 
pairs, exact boundaries of quantum phase transitions 
among various competing states and the number of the 
collective modes. In the regime where accurate Monte 
Carlo simulations can be carried out without the sign 
problem, detailed quantitative comparisons with exper- 
iments are possible, including the quantum phase tran- 
sition from the AF to the SC phases as a function of 
doping. 
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